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Improved Method for the Timed Dosage of Medicaments 

The invention relates to a method for the controlled dosage of a medicament as a function of time 
by means of a method for the determination of a corresponding dosage profile and corresponding 
control of a dosage instrument. 

5 Computer-controlled infusion pumps, the delivery functions of which are determined by means of 
a pharmacokinetic model, are known according to the prior art by the keyword "TCI" (= Target 
Controlled Infusion) and are commercially available. The main application field of TCI is the 
control of intravenously administered narcotics (for example propofol, marketed as Diprifusor™ 
by AstraZeneca (Product information "Diprifusor™: Target Controlled Infusion (TCI) in 

10 anaesthetic practice", AstraZeneca Anaesthesia, New Edition (1998)). A disadvantage of these 
known methods is that the pharmacokinetic model is a three-compartment model fitted to 
experimental plasma data. With such a "black-box" method, there is no opportunity for the 
patient's individual physiological factors to be taken into account in the pharmacokinetic model. In 
contrast to this, physiology-based PK/PD models such as PK-Sim® developed by Bayer 

15 Technology Services GmbH (www.PK-Sim.com ; S. Willmann, J. Lippert, M. Sevestre, J. 
Solodenko, F. Fois, W. Schmitt: "PK-Sim®: a physiologically based pharmacokinetic 'whole -body' 
model", Biosilico 1, 121-124 (2003)), makes it possible to describe the influence of individual 
physiological and anatomical parameters such as organ size and composition, blood flow rates, etc. 
on the pharmacokinetic behavior of medicaments as a function of time (for example DE App. No. 

20 10160270 and DE App. No. 10345836). These physiological and anatomical parameters can in 
turn be attributed to a few easily measurable quantities such as body weight and body mass index. 
DE App. No. 10345837 furthermore describes the way in which biochemical and genetic 
information, for example expression data of metabolically active enzymes or active transporters, 
can also be employed in order to determine a dose individually adapted to the patient. These 

25 systems are pure simulation models, however, which do not allow direct support of an application 
to the patient or clinical applications. 

On the basis of said prior art, it is an object of the invention to develop an improved method which 
permits exact dosage of a medicament as a function of time while taking into account the patient's 
individual physiological, anatomical, biochemical and genetic factors. 

30 The exact dosage of a medicament as a function of time is crucial for the safety and success of the 
treatment in many indication fields (for example anesthesia, diabetes, shock, sepsis, cardiovascular 
failure, asthma, cancer). With the aid of electronically controlled infusion pumps, medicaments 
can be administered with an arbitrarily predetermined time -variable rate. The resulting 
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concentration-time profile and effect-time profile do not depend only on the dosage profile, 
however, but are essentially determined by the pharmacokinetic (PK) and pharmacodynamic (PD) 
properties of the substance in question. Physiology-based pharmacokinetic (PBPK) and 
pharmacodynamic computer models are in principle capable of simulating the concentration-time 
5 profile as well as the effect-time profile of a chemical substance in a patient's body, while taking 
into account individual physiological, anatomical, biochemical and genetic parameters. 

The invention therefore relates to a method in which a PBPK/PD model is used in order to 
determine an optimal time profile of the dose for the individual patient by iterative adaptation 
either of the concentration-time profile, in the plasma or at the target site, or of the 
10 pharmacodynamic effect-time profile to a predetermined target time profile. The optimized time 
profile of the administration is then used as the input function for a dosage device. In combination 
with real-time measurements of physiological parameters, it is possible to set up a closed control 
loop that is far superior to the black-box methods known in the prior art, which do not utilize 
physiological knowledge. 

15 It is particularly advantageous to combine the method with a real-time measurement of 
physiological parameters which may be subject to variations as a function of time, as is the case 
for example within the expression of metabolically active enzymes. 

The essential feature of the invention consists in the combination of a physiology-based 
pharmacokinetic and/or pharmacodynamic (PBPK/PD) model with an automated dosage device, 
20 for example an electronically controlled infusion pump. PBPK/PD models are advantageous over 
non-physiological compartment models since they are capable of describing the effect of 
individual physiological, anatomical, biochemical and genetic factors on the pharmacokinetics and 
-dynamics in more detail. 

By using the model PBK/PD, it is surprisingly now possible to determine very individual dosage 
25 rates adapted to the patient. 

A schematic representation of the method according to the invention is shown in Figure 1. The 
main component is a PBPK/PD model (101) which describes a particular type of body, particularly 
preferably that of a mammal, and requires a series of different parameters as input quantities: 

1.) Substance-specific parameters of the medicament to be administered (102). Typical 
30 substance-specific parameters are e.g. physicochemical parameters such as lipophilicity, 

binding constants to human serum albumin and/or other plasma proteins, unbound plasma 
fraction, solubility in aqueous buffer solution or in intestinal fluid, size of the molecule 
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(expressed by the molecular weight or molar volume), hepatic and/or renal clearance, 
permeability coefficients e.g. across artificial or biological membranes, and equilibrium 
distribution coefficients between plasma (or blood) and the various organs. 

2. ) Species-specific physiological, anatomical, biochemical and/or genetic input parameters, 

which are characteristic of the patient in question (103). These types of parameters include 
in particular body weight, volume components of individual organs with respect to the 
entire body volume, blood flow rates of individual organs, water, fat and lipid components 
of the individual organs, and parameters which characterize the expression and function of 
metabolically active enzymes (particularly in the liver and intestines) or the expression and 
function of proteins for the active transport of molecules through cell membranes. These 
parameters either may be directly measured or are correlated for particular populations 
with easily determinable patient parameters such as age, sex, body weight and lean body 
mass (P. S. Price, R. B. Conolly, C. F. Chaisson, E. A. Gross, J. S. Young, E. T. Mathis, D. 
R. Tedder: "Modeling interindividual variation in physiological factors used in PBPK 
models of humans", Crit. Rev. Toxicol 33, 469-503 (2003)). Individual parameters or 
several parameters among these may be used in the method according to the invention. 

3. ) A dosage profile which describes the dose administered as a function of time (104). 
The following is also specified: 

4. ) A target profile which indicates the intended concentration-time profile of the medicament 

to be administered in plasma, blood or directly at the biochemical target in the target organ 
or the intended effect-time profile ("SETPOINT profile", 105). This target profile is 
indication-specific and dependent on the medicament in question. In many clinically 
relevant cases, for example the administration of narcotics, a rectangular effect profile 
with maximally steep edges is desired, i.e. the intended effect (here narcosis) is to set in 
spontaneously then remain as constant as possible over a defined time period, and then 
drop off again rapidly at the end of the treatment. Either the time profile may be a simple 
function of time Z(t) or a tolerance range may alternatively or additionally be specified 
(defined as an interval between a maximum value and a minimum value [Z min (t)..Z max (t)]). 

On the basis of a start function to be selected expediently for the dosage profile, the PBPK/PD 
model calculates the individual concentration-time profile or effect-time profile for the substance 
in question with the aid of the information from 1.) - 2.) ("ACTUAL profile", 106). 
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Here, fundamentally different routes for the administration of the medicament may be envisaged in 
the PBPK/PD model. The medicament may be administered intravenously, intra-arterially, 
intraperitoneally, intramuscularly, subcutaneously, topically, orally or inhalatively via the upper 
airways and/or the lungs. Intravenous application is particularly important and preferred. 

The next step of the method consists of an iterative optimization process, in which the dosage 
profile is varied until the simulated concentration-time profile or effect-time profile matches the 
SETPOINT profile (107, 108). As a result of this numerical optimization, the dosage time profile 
is obtained which causes the desired concentration-time or effect-time profile for the relevant 
substance in the individual patient, or exhibits the least deviation therefrom. Gradient methods, 
gradient-free methods, stochastic methods or evolutionary methods may be envisaged as numerical 
optimization methods. Among gradient methods, the quasi-Newton or Newton methods are 
particularly preferred, and in particular the nested interval method among gradient-free methods. 
The Monte-Carlo method is preferred in particular among stochastic methods, and the method of 
genetic optimization represents a particularly preferred form of an evolutionary method. 

This dosage profile is employed in the last step in order to control an automatic dosage device. 

Accordingly, the present invention also relates to a method for the controlled dosage of a 
medicament as a function of time, consisting of the following steps: 

a) Specification of an indication- and substance-dependent target profile, which 
indicates a desired concentration-time profile or a desired effect-time profile, 

b) Physiology-based pharmacokinetic and/or pharmacodynamic simulation with a time- 
variable application profile while taking into account individual anatomical, 
physiological and/or genetic parameters of the body to be treated and substance- 
specific input parameters of the medicament to be administered, 

c) Iterative numerical adaptation of the application profile until the simulated time 
profile matches the predetermined target profile, or shows the maximum achievable 
match therewith, 

d) Control of a dosage device on the basis of the result in c). 

The method according to the invention is particularly advantageous when physiological 
parameters, which have an effect on the pharmacokinetic or pharmacodynamic behavior of the 
medicaments, are time -variable. Examples of physiological parameters which can change during 
the therapy method are the blood flows of the gastrointestinal tract and the liver, which are 
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influenced and modified inter alia by the absorption of nutrients, or (over a longer timescale) the 
expression rates of transporters or metabolically active enzymes. If such parameters are measured 
in a real-time during the therapy, then their time-varying effect can also be taken into account in 
the optimization of the dosage profile by the method according to the invention. 

5 In certain cases, for example anesthesia, it may be expedient to not only regulate the dosage device 
by the desired pharmacokinetic or pharmacodynamic target profile, but also to employ an external 
measurement quantity as a further input parameter. 

In a particular embodiment, the success of the therapy is additionally monitored online by one or 
more suitable measurement probes, which measure for example the depth of narcosis in anesthesia, 

10 and the measurement signals are integrated as additional input quantities in the method. The 
control of the dosage device is then not purely regulated by the pharmacodynamic or 
pharmacokinetic target profile, but also incorporates external measurement signals. For example, 
the supply of a narcotic may be increased when the measured depth of narcosis falls below a 
critical value. In this way, it is possible to set up a closed control loop which optimizes the time 

15 dosage of the medicament by incorporating online measurement values and physiological 
simulations in real-time. In this embodiment, the response of the pharmacodynamic or 
pharmacokinetic profile to changes in the application rate, which is known from the physiological 
simulation, is used in order to readapt the rate of the medicament supply to the patient's 
requirements as indicated by the measurement probe. In another particular embodiment, only the 

20 measurement probe signal is employed temporarily for the regulation if there is a critical situation 
(for example if the patient is about to wake early from the narcosis). 

All methods based on said parameters are in principle suitable as simulation methods, the methods 
claimed in DE App. No. 10160270 and DE App. No. 10345836 being particularly suitable and 
preferred according to the invention. 

25 Besides the application of the method according to the invention as an aid for carrying out a 
medical therapy, the method according to the invention may also be used as an aid directly in 
clinical trials or animal trials, for example in order to start off the runs with clinically "expedient" 
doses and to minimize the typical "settling in" of the doses, i.e. the empirical-iterative arrival at 
excessive or insufficient doses which alternatingly approach the optimum, and therefore minimize 

30 the burden on the bodies being treated and maximize the likelihood of the experiment's success. 



Humans and animals are therefore suitable as a target group for the application of the method 
according to the invention, especially humans and economically useful, breeding, laboratory, test 
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and pet animals. The method is more particularly preferably used as an aid for the therapeutic 
treatment of humans or clinical trials on humans. 

Economically useful and breeding animals include mammals, for example cows, horses, sheep, 
pigs, goats, camels, water buffalo, donkeys, rabbits, fallow deer, reindeer, animals prized for fur, 
5 for example mink, chinchillas, raccoons, birds, for example chickens, geese, turkeys, ducks, 
pigeons, and bird species kept at home and in zoos. 

Laboratory and test animals include mice, rats, guinea pigs, hamsters, rabbits, dogs, cats, pigs and 
apes, respectively in all species, subspecies and breeds. 

Pet animals include dogs and cats. 

10 Electronically controlled infusion pumps, inhalers or electronically controlled release capsules for 
oral application may be envisaged as the electronically controlled dosage device. Electronically 
controlled infusion pumps are particularly suitable. 

The method according to the invention is particularly advantageous for those indications and 
active agents which have only a narrow "therapeutic window". A narrow therapeutic window 

15 means that there is only a small concentration range in which the desired pharmacological effects 
do actually occur but at the same time no undesired side-effects are to be observed. Examples of 
indication fields with narrow therapeutic windows are all types of cancer diseases, infectious 
diseases, in particular bacterial and viral infections, cardiovascular diseases, in particular high 
blood pressure, lipidemia, angina pectoris and myocardial infarction, diseases of the central 

20 nervous system such as Alzheimer's disease, schizophrenia, epilepsy, chronic headaches 
(migraines), analgesia and anesthesia, psychiatric diseases, in particular depression and anxiety, 
metabolic diseases, for example diabetes and impairments of fat metabolism (obesity), respiratory 
diseases such as asthma and bronchitis, immune diseases, in particular allergies, rheumatism and 
multiple sclerosis, diseases of the gastrointestinal tract, in particular ulcers of the stomach and 

25 duodenum and Crohn's disease, as well as vascular diseases, in particular those which cause 
erectile dysfunction, and states of acute shock. 
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Example 

The invention will be explained in more detail below with the aid of an example; this example 
illustrates the present invention but is in no way to be regarded as a limitation. 

The example is based on simulations with the physiology-based pharmacokinetic model PK-Sim® 
5 developed by Bayer Technology Services. 

Theophylline, a medicament for asthma therapy, was selected as an example substance. The 
peculiarity of this substance is that it has only a very narrow concentration range in which the 
substance is therapeutically effective. This so-called "therapeutic window" is described in the 
literature as being 5-10 to 20 mg/1 (plasma concentration) (A. Ohnishi: "A Review of Clinical Use 

10 of Theophylline in Acute Asthma: Factors Influencing Kinetic Disposition and Drug Interactions", 
Methods Find. Exp. Clin. Pharmacol. 22, 253-258 (2000)). A sufficient anti-asthmatic effect is not 
generally achieved at plasma concentrations below 5 mg/1, and the risk of undesired side-effects 
such as nausea, trembling and accelerated heart rate increases above a plasma concentration of 20 
mg/1 (A. Ohnishi: "A Review of Clinical Use of Theophylline in Acute Asthma: Factors 

15 Influencing Kinetic Disposition and Drug Interactions", Methods Find. Exp. Clin. Pharmacol. 22, 
253-258 (2000)). An added difficulty in the therapeutic use of theophylline is the known fact that 
the pharmacokinetics of this medicament depend strongly on the patient's individual factors, for 
example age, sex or body fat level. Particularly the metabolization rate (P. A. Mitenko, R. I. 
Ogilvie: "Pharmacokinetics of intravenous Theophylline", Clin. Pharm. Therap. 14, 509-513 

20 (1973); J. P. Jameson, A. Munyika: "Theophylline Pharmacokinetics in Black Zimbabwean 
Males", Therap. Drug Monitoring 12, 54-58 (1990); P. Gal, W. J. Jusko, A. M. Yurchak, B. A. 
Franklin: "Theophylline disposition in obesity", Clin. Pharmacol. Ther. 23, 438-44 (1978); S. H. 
D. Jackson, A. Johnston, R. Woollard, P. Turner: "The Relationship between Theophylline 
Clearance and Age in Adult Life", Eur. J. Clin. Pharmacol. 36, 29-34 (1989)) and the free plasma 

25 fraction (Z. Trnavska: "Theophylline Protein Binding", Arzneim.-Forsch. 40, 166-170 (1990); Z. 
Kato, O. Fukutomi, N. Kondo: "Developmental Changes of unbound Theophylline", Ann Allergy 
Asthma Immunol. 80, 517 (1998)) are individually different and consequently lead to individually 
different plasma concentration profiles. In clinical use, this problem is tackled by experimentally 
monitoring the plasma concentration of theophylline (A. Ohnishi: "A Review of Clinical Use of 

30 Theophylline in Acute Asthma: Factors Influencing Kinetic Disposition and Drug Interactions", 
Methods Find. Exp. Clin. Pharmacol. 22, 253-258 (2000)). This method is comparatively 
elaborate, it requires repeated blood samples to be taken from the patient and is furthermore only 
expediently usable for long-term therapy with theophylline. 
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The way in which the problem of an optimal theophylline dose tailored individually to the patient 
can be resolved by the method according to the invention will be described below: 

The following values for theophylline are used as substance-dependent input parameters (Table 1): 

Table 1 



PARAMETER VALUE 



Lipophilicity (LogMA) 


-0.2 


Free plasma fraction 


59% 


Solubility 


2800 mg/1 


Molar mass 


180 daltons 


Intestinal permeability 


2.97 * 10" 6 cm/s 


(only for taking orally, Fig. 2c) 





5 The value of 59% used for the free plasma fraction represents an average value for an adult human 
(Z. Trnavska: "Theophylline Protein Binding", Arzneim.-Forsch. 40, 166-170 (1990); Z. Kato, O. 
Fukutomi, N. Kondo: "Developmental Changes of unbound Theophylline", Ann Allergy Asthma 
Immunol. 80, 517 (1998)). The following physiological and anatomical parameters were also 
assumed for the calculations with PK-Sim® (default parameters for humans, Table 2): 

10 Table 2 



Organ/compartment 


Volumes [ml] 


Blood 


Venous blood pool 


250 


-/- 


Arterial blood pool 


140 


-/- 


Lungs 


670 


4670 


Stomach 


150 


60 


Small intestine 


640 


600 


Large intestine 


370 


240 


Pancreas 


100 


60 


Spleen 


180 


180 


Liver 


1710 


390 


Gall bladder 


20 


0 


Kidneys 


720 


1133 



Blood flow rates [ml/min] 
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Brain 


1 A 


7AA 


Heart 




Z<4U 


Muscle 






Bone 


12060 


167 


Skin 


3020 


50 


Fat 


10060 


300 


Testicles 


35 


2.6 



It will first be shown that in principle it is possible to describe the plasma concentration-time 
profile of theophylline by PK-Sim® correctly on the basis of the aforementioned values. To this 
end, simulation results were compared with experimentally measured plasma concentration-time 
profiles (S. Valenti, P. Crimi, V. Brusasco: "Bioavailability and Pharmacokinetics of a New 
5 Controlled Release Theophylline Preparation in the Form of Capsules Containing Pellets", Respi- 
ration 52, 195-200 (1987); M.Muller, B. v. Osten, R. Schmid, E. Piegler, I. Gerngross, H. Buch- 
egger, H. G. Eichler: "Theophylline kinetics in peripheral tissues in vivo in humans", Naunyn- 
Schmiedeberg's Arch. Pharmacol. 352, 438-441 (1995); M. C. Meyer, E. J. Jarvi, A. B. Straughn, 

F. R. Pelsor, R. L. Williams, V. P. Shah: "Bioequivalence of Immediate -release Theophylline 
10 Capsules", Biopharm. Drug Dispos. 20, 417-419 (1999)) and a muscle concentration-time profile 

from the literature (M.Muller, B. v. Osten, R. Schmid, E. Piegler, I. Gerngross, H. Buchegger, H. 

G. Eichler: "Theophylline kinetics in peripheral tissues in vivo in humans", Naunyn- 
Schmiedeberg's Arch. Pharmacol 352, 438-441 (1995)) (Figure 2). As can be seen in Figure 2, the 
match of the simulated curves with the experimentally determined curves is very good. 

15 Two virtual individuals will be considered below for the optimization of individual dosage 
profiles, who differ by way of example in respect of their total metabolization rate ("clearance"). 
Individual A has a total plasma clearance (CL ) of 0.4 ml/min/kg, and individual B has 
1.0 ml/min/kg. These values lie within the typical variation range of the plasma clearance of 
theophylline in humans (P. A. Mitenko, R. I. Ogilvie: "Pharmacokinetics of intravenous 

20 Theophylline", Clin. Pharm. Therap. 14, 509-513 (1973); J. P. Jameson, A. Munyika: 
"Theophylline Pharmacokinetics in Black Zimbabwean Males", Therap. Drug Monitoring 12, 54- 
58 (1990); P. Gal, W. J. Jusko, A. M. Yurchak, B. A. Franklin: " Theophylline disposition in 
obesity", Clin. Pharmacol. Ther. 23, 438-44 (1978); S. H. D. Jackson, A. Johnston, R. Woollard, P. 
Turner: "The Relationship between Theophylline Clearance and Age in Adult Life", Eur. J. Clin. 

25 Pharmacol. 36, 29-34 (1989)). The other relevant physiological and anatomical parameters, which 
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are used in the following simulations, are assumed to be equal for both individuals and identical to 
the values collated in Table 2. At this point, it is naturally also conceivable to take into account a 
multiplicity of other individual differences, for example in free plasma fraction, body weight, fat 
level body or blood flow rate in one or more organs. 

A plasma concentration of 15 mg/1 is selected as the SETPOINT profile in this example, which is 
intended to be kept constant over a time period of 12 hours. The selected tolerance range is 
± 1 mg/1 (represented in Figure 3), i.e. 



2^(0 = 



Z(t) = 



0 , t<0 
15mg/L , 0<t<12h 
0 , t>12h 



0 , t<0 
14mg/L , 0<t<12h 
0 , t>12h 



and Zmax (0 = 



and 

0 , t<0 
16mg/L , 0<t<12h 
0 , t>12h 



(1) 



Figure 4 shows the plasma concentration-time profiles simulated by PK-Sim® for the two 
individuals after a 1 mg/kg intravenous bolus administration of theophylline. The curves are used 
10 as a starting point for optimization of the dosage profile. 

The infusion rate, i.e. the intravenously applied dose per time interval, is then varied iteratively in 
the next step. The increment of the dosage interval in this example is 0.1 hour = 6 minutes. This 
value is expediently to be selected so that it is adapted to the distribution and elimination kinetics 
of the substance. 



15 The simulation result obtained (= ACTUAL profile) is then compared with the SETPOINT profile. 
In the case of a deviation of the ACTUAL profile from the SETPOINT profile which lies outside 
the valid tolerance range, the amount applied in the corresponding time step is varied until the 
ACTUAL profile matches the SETPOINT profile overall in the scope of the tolerance range. The 
optimization of the application profile was carried out by manual iteration in the present example. 

20 In the scope of a genuine therapy optimization, suitable numerical optimization methods are to be 
used at this point. Gradient methods, in particular quasi-Newton or Newton methods, as well as 
gradient-free methods such as nested intervals and stochastic methods such as Monte-Carlo 
methods may for example be envisaged here as suitable optimization methods. 

A respective result of such a manual iterative optimization process is represented in Figures 5 (a) 
25 and (b). The dosage profiles associated with these ACTUAL profiles are shown in Figures 6 (a) 
and (b). These dosage profiles determined in this way are used in the last step of the method as an 
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input function for the control of a conventional dosage machine, for example an electronically 
controlled infusion pump. 
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Description of the Figures 

Figure 1: Schematic representation of the method for the timed dosage medicaments. 

Figure 2: Comparison between experimental concentration-time profiles of theophylline from the 
literature with simulation results by PK-Sim® with metabolization rates adapted to the respective 
5 experimental data: 

a) Intravenous infusion of 5.58 mg/kg over 20 minutes. The plasma clearance here is 
CL = 0.40 ml/min/kg. (source: [12]) 

b) Intravenous infusion of 240 mg over 10 minutes. The plasma clearance here is 
CL = 0.25 ml/min/kg. (source: [13]) 

10 c) Peroral administration of 200 mg of an immediately releasing formulation (mean and 
standard deviation of three bioequivalent immediate release capsules [14]). The plasma 
clearance here is CL = 0.50 ml/min/kg. 

Figure 3: Specification of the target profile (= SETPOINT profile) 

Figure 4: Simulation result for the intravenous bolus administration of 1 mg/kg theophylline in 
15 individuals A (CL = 0.4 ml/min/kg) and B (CL = 1.0 ml/min/kg). 

Figure 5: Simulated plasma concentration-time profile (ACTUAL profile) after iterative 
adaptation of the dosage profile for individuals A (a) and B (b). 

Figure 6: Optimized dosage profile for individuals A (a) and B (b). 



